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Abstract. Motivated by a recent use of Glauber dynamics for Monte-Carlo simulations of path 
integral representation of quantum spin models [4] , we analyse a natural Glauber dynamics for 
the quantum Ising model with a transverse field on a finite graph G. We establish strict mono- 
tonicity properties of the equilibrium distribution and we extend (and improve) the censoring 
inequality of Peres and Winkler to the quantum setting. Then we consider the case when G is 
a regular fe-ary tree and prove the same fast mixing results established in [7] for the classical 
Ising model. Our main tool is an inductive relation between conditional marginals (known as 
the "cavity equation") together with sharp bounds on the operator norm of the derivative at 
the stable fixed point. It is here that the main difference between the quantum and the classical 
case appear, as the cavity equation is formulated here in an infinite dimensional vector space, 
whereas in the classical case marginals belong to a one-dimensional space. 
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1. Introduction 

In the last years there has been a somewhat intense research in probabihstic representation of 
quantum spin systems with very interesting progresses exploiting stochastic geometry methods. 
We refer the interested reader to [2j and references therein. More recently, physicists started to 
consider quantum Heat Bath algorithms applied to the imaginary time path integral representa- 
tion of e.g. the quantum ising model in a transverse field (see [5] and references therein). From 
a mathematical point of view that accounts to introduce a so-called Glauber dynamics for the 
dynamical variables appearing in the path integral representation of the model, a subject that, 
in the context of classical spin models, has been thoroughly investigated in the mathematical 
literature in the past fifteen years because of its many different facets and connections with dif- 
ferent areas in mathematics, physics and theoretical computer science. In the quantum setting 
very little is known rigorously and we are only aware of a recent contribution [IJ where, for a 
particular choice of the Glauber chain, exponentially fast approach to the imaginary time equi- 
librium measure is proved via a kind of Bochner-Bakry-Emery method in a high temperature 
regime. 

Here we investigate the simple heat bath dynamics for the imaginary time path integral 
representation of the quantum Ising model with a transverse field on a finite connected graph G. 
We first focus on properties of the single site measure and establish key monotonicity properties 
as a function of the longitudinal magnetic field and of the boundary conditions. A remarkable 
consequence of monotonicity is the validity of the so-called Peres- Winkler censoring inequality 
[8] which we re-establish and extend to the quantum case. Then, motivated by the analysis of 
the same model made by physicists on the Bethe lattice [4], we specialize to the case when the 
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graph G is a regular b-ary tree and we succeed in extending to the quantum case much of the 
analysis made in [7J for the classical setting. More precisely we prove that, for any boundary 
condition r, there exists exponents k = k{t), 7 < 1 such that, if Kjb < 1, then the Glauber 
chain on a finite subtree of depth i with boundary condtion r at the leaves, has a spectral 
gap uniformly positive in i. In the same uniformity was proved for the logarithmic Sobolev 
constant and, as a consequence, the mixing time cound not grow faster than linear in £. In the 
quantum setting the logarithmic Sobolev constant is easily seen to be infinite. Nevertheless, 
using our extension of censoring and always assuming K76 < 1, we prove that the mixing time 
grows linearly in i. Finally, in complete analogy with the fast mixing result inside the pure plus 
phase for the classical Ising model, we show that, if the boundary condition is the homogeneous 
plus configuration then the key exponent k satisfies k ^ 1/6 for all values of the thermodynamic 
parameters. Therefore, in this case, K76 < 1. The same holds for all boundary conditions when 
the thermodynamic parameters are in the uniqueness regime. 

All the main technical efforts in the paper are concentrated on bounding the exponent k and 
it is at this point that the analysis of the quantum case becomes much more involved than its 
classical analog. Computing the exponent k requires in fact a detailed study of the quantum 
cavity equation for the one site conditional marginal distributions and, in particular, of the 
derivative operator at the fixed point corresponding to the "plus" phase. In the classical case 
such marginals are parametrized by just one number and therefore the corresponding cavity 
equation is just a recursive non-linear equation in R. In the quantum case instead the cavity 
equation becomes a recursive functional equation in a infinite dimensional vector space. 
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2. The quantum Ising model 



2.1. The matrix formulation. We consider a finite graph G = {V,E) with vertex set V and 
edge set E. A classical spin configuration on y is o" = {ai)i^v S Now we define Ti, a 
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2l^l-dimensional Hilbert space, by giving the orthonormal basis {\a))^^^j_^'^v ■ We also denote 
{a\ the transpose of any basis vector \a). Given i V, we define two linear operators on Ti by 
their action on the orthonormal basis: 

(T||a) = (Tiler) 

<k) = 

where cr* is the same spin configuration as a, except for the spin at i which is flipped. Given a 
set of parameters (3,h, X (/3 > is the inverse temperature, /i G R the longitudinal and A ^ 
the transverse field) we define the Hamiltonian operator 

e={i,j}eE iev ieV 

and define the average of an observable O (a linear operator on H) at inverse temperature (3 as 



Tr (O 



Tr(e- 



2.2. The path-integral representation of the quantum Ising model. We introduce first 
a reference measure <pi^x that corresponds to a single free spin on some interval / C R. A single 
spin is a cadlag function <t : / — )• {±1}. We call S/ the set of single spins, and endow it with 
the usual Skorohod topology. Note that two spin configurations are close when the have the 
same initial value, same number of flips and close time location for the flips. This set is a Polish 
space as the set of spin configurations with jumps at fractional times is countable and dense in 
S. The corresponding cx-algebra is generated by the events {cr (t) = +}, for t £ I. We will also 
consider the L^-distance 



\cT-a'\ 



j\a{t)-a' {t)\dt. 



between two spin configurations. Now we define the distribution (pi^x as follows: we consider a 
Poisson point process with intensity A on /, and define <t, with equal probability, as one of the 
two configurations that flip exactly at the points of the Poisson process. 

Now we introduce spin interactions and external fields. Given a bounded interval / and two 
bounded fonctions cr, r : / ^ R defined Lebesgue a.s., we call 

<T»T = j a{t)T{t)dt. (2.1) 

When I is clear from the context (usually / = [0,/3]), we drop it from the notation. Note that 
the product defined by (12. ip is additive in the sense that, for /i,/2 disjoint, 

a • T = a • T + a • T . 

huh h h 

A spin configuration is <t : y x / — )• {±1} such that CTj G S/, for all i £ V. Consider some 
/? > (inverse temperature), /i : F x [0, /3] — > R integrable (the longitudinal field), A ^ (the 
transverse field). We define two probability measures on the set of spins configurations: 



fJ'G;P,h,x{da) = exp j ^ CTj • cr^ + ^ /ij • cTj j J| (^[o,/3],A(dcri) 




n '/'[0,/3],A(dcri 
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where ^G;/3,h.,A and Z^^^ ^ ^ are the appropriate normahzation constant that turn fJ,G;i3,h,x and 
/^G-^ h A ^^^^ probabihty measures on the set of spin configurations. They exists as the argument 
of the exponential is triviahy bounded by 13 \E\ + /3 \\h\\-^ \V\. The boundary condition cr (0) = 
a {13) is cahed periodic imaginary time boundary condition. Note that we could also allow 
A : y X / ^ R+. 

The above measure are often approximated, in the physics literature, by a discrete Ising model 
with vertex set F x {0, . . . , N} and edge set E x {0, . . . , N}UV x {(0, 1) , . . . , {N - 1, N) , {N, 0)}. 
This is usually referred to as a consequence of the Suzuki- Trotter formula. The strength of 
interaction is different for transversal edges than for longitudinal edges. Although we do not use 
this approximation (which corresponds to enabling the discontinuities only at times /3k /N), it 
explains that the above two measures share most of the properties of discrete Ising models. 

Proposition 2.1. [DLR Equation]. Let I = [/3i,/32] C [0, /3] be an interval and A C V. Let 
^ f^G:i3,h,x- Conditionally on cr^(^y^y.jy , the distribution of cr^y^^j is the quantum Ising model 
l^{A,Ey,i,h',\ {-Wa (/3i) ,<ta {132)) with external field h' = h + J2ieA,j^A^j^i- 

Proof. This is a consequence of the independence of the restrictions of the Poisson point process 
to different regions together with the additivity of the product (12. ip . □ 



Proposition 2.2. Let /3 > 0, /i € R and A ^ 0. Then, for any observable O, 

inx /iy;AM((^(/3)|Ok(0))) „^ 

while, for any diagonal observable O, 

(0)/3,M = ^^V^p^H^{W)\o\cT{m■ (2.3) 

Proof. We define a matrix Wp^h,\ by the prescription of its coordinates 



= 2'^' j exp j ^o-i • cTj + 1 • cTj j l|^(o)=<^,^(^)=^/} JJ dv9[o,/3],A(o-i)- 
Once we prove that 

Wp,h,x = e"'^^'^' X exp {-f3H) (2.4) 
the claim follows from the definition of the average of an observable: 

\^)is,h,x Xr(e-'3^^) EM^'^^W) 
_ Ea,a'{^\OW'){cr'\e-^''W) 

as 

{a'\e-^^\a) = e^'\''\{a'\Wp,H,,\a) 
= e^^\^\2\^\Zv;px>^ //y;/3,/.,A(T(0) = a',a{P) = a). 

So we focus now on the proof of (j2.4p . According to the additivity of the product (12. ip and to 
the independence of the restriction of the Poisson process to [0,/3i] and [/3i,P2], we have 

Wp,+p,,h,x = Wft,MxT^;32,M,V/3i,/32^0. (2.5) 

Now we compute the asymptotic of {o'\Wi3^h,x\(^') when (3^0. The probability that the Poisson 
point process has two points or more on y x [0, f3] is 0(/3^), therefore when a and a' differ at two 
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points or more, {o'\Wis^h,x\(^') = 0{(3'^). When a' = cr*, the probabihty that cTj is discontinuous 
is /3A + 0(/32), therefore {a\W[s,h,xW') = PX + 0{^^) = p{a\H\a') + 0(/32). Finahy, when 
o"' = (T, up to 0(/3^) only the configurations with no flips (probability exp(— /3A|y|)) contribute 
to {cr\Wp^h,\W) , therefore in that case 

{(r\Wp^h,xW') = exp I /3 I ^aiCT,- + J^a^ | - p\\v} 



= I- P{a\H\(T') - pX\V\ + 0{p'^). 

In other words we have shown that 

Wp,h,x = {l-p\\V\)I-PH + 0{p^) 

as /3 — >■ 0, where / is the matrix that corresponds to the identity of T-L. Combining ()2.5p and 
(j2.6p we obtain, for any /3 ^ and any cr', 



(2.6) 



dp' 



/3'=0 



= {-X\V\I-H)xW^,h,xW'). 
This proves (j2.4p as exp(— /3A|y |/ — /3ff) is the unique solution to this differential equation. □ 



3. The single site spin measure 

In all this section we assume that V = {x} is a single site. Most of the time /3, A and sometimes 
h will be from the context and they will not appear in the notation. We state several important 
facts about this single site measure. In what follows || • ||tv denotes the usual variation distance 
(see also (fO]) . 

Proposition 3.1. Let f3,M < oo. There exists 7 < 1 and T = l/log(3) such that, for any 
h,h',X such that \\h\\i, \\h'\\i,X < M, for any vr S {0,per}, the following inequalities holds 



I^;^'-m^IItv ^ m'-hw,. 



(3.1) 
(3.2) 



Proof. Inequality (|3.ip is an immediate consequence of the fact that, under our assumptions, 
/^h,({+}) is uniformly bounded from below. For (|3.2p we introduce cp'^ , the single spin free 
measure conditioned on the imaginary time boundary condition vr and compute 



\l^h' ^^11 TV 



< - 



1/ 

2 



/e^'"^'d99^(cr') /e^"^'dv9^(cr') 

dv9^(cr)dv9^(cr' 



»cr+h'mcr> 



dip'^ {cr)dip'^ {a') 



2 



h»cr+h •a' 



,{h'-h),(cT-cr') _ ^ 



dv9'^(cr)d^^(cr') 



•(T+h'»a' 



d^p-^ {a)dip'^ {a') 



exp(||^' - h\\i) - 1 



The claim follows from the inequality (e^ — 1) /2 ^ Fx, for every x > such that Fx ^1. □ 
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The second property is strict monotonicity , a crucial fact for our analysis of the dynamics. 
Let / : S — )• R be a function of a single spin configuration. We say that / is increasing 
when cr (t) ^ cr' (t) for almost every t e [0,/3] implies / (cr) ^ / (c')- When are two 
probability measures on S such that, for all increasing function /, i' (/) ^ i^' (/), we say that i' 
is stochastically smaller than u' and we will write v u' . 

Theorem 3.2. Let I3,M < oo. There exists c > such that, for any h,h',X such that 
A < M, for any vr G {0,per}, if h' ^ h point-wise a.s. then fij^, >~ fij^. More- 
over, for any f increasing, 

f^l'in-f^lif) > c{f{+)-f{-)){h'-h).l. (3.3) 

A consequence of monotonicity is the FKG inequality: 

Corollary 3.3. Let /3,A < oo and h with ||^||oo < oo. Then, for any it G {0, per}, any f,g 
increasing, 

f^Hfa) ^ (3.4) 

The main idea in the proof of Theorem 13.21 is an explicit coupling of two single site measures 
with the same parameters but different imaginary time boundary condition. 

Lemma 3.4. Fix /3,A > and h G L^([0,/3]). Consider two independent spins variables 
(T^ ~ /i(.|<T(/3) = +) and ~ |u(.|<t(/3) = — ). Call T^ (respectively T~ ) the last flipping time 
(0 if none) of (respectively of cr^}, and T = max (T+, T"). Consider the joint distribution 
^ on (<T+,cr^) as follows: 

(1) Let = + and ""(^^j = — • 

(2) Take cjjj^^ = ''■[q^) ^ l^i-Wi"^) = ^) where e = + if T^ < T~ , e = — otherwise. 
Then ^ is a monotone coupling of fi{.\cr{l3) = +) and fj,{.\a{f3) = —). It satisfies 

$(cr+ = +, cr- = -) = ^(cr = +[ct(/3) = +) x /i(cr = -\a{p) = -). (3.5) 



Proof of Lemma \3.4\ ^ is monotone by construction (cr+ ^ a point-wise a.s.) while (|3.5p 
is a consequence of cr^ = ± <^4> T^ = T^ = 0. In order to prove that ^ is a coupling we 
need to check that it has the correct marginals. Let t G (0, /?). Conditioning on T^ < t is 
the same as conditioning the Poisson point process on having no point in (t, /?) and therefore 
a ~ /i(.|cr(/3) = +,r+ < t), restricted to (0,i), has distribution ^^^^ ;^(.|<T(i) = +). Conditioning 
on = t is more delicate as this is a zero probability event, still it has a very precise meaning 
in terms of the Poisson point process, since it requires that there is a point at t and not point in 
{t,l3). It is well known that the conditional distribution of the points in {0,t) is an independent 
Poisson point process with the same intensity. Therefore, cr ~ fi{.\cr{(3) = +,T+ = t), restricted 
to (0,t), has distribution fj.t,h,x{-\f^{t) = —)■ The same holds if we replace + by — and vice-versa. 
Since T~^ and T~ are independent, this proves that \I' is indeed a coupling of //(.|<t(/3) = +) and 
M(-k(/3) = -). □ 

Proof of Theorem \3.S\ We begin with the case without periodic imaginary time boundary con- 
ditions, that is vr = 0. It is enough to quantify the influence of a unitary increase of the field, so 
we consider A/i : [0,/?] R"*". We remark that 



-^fJ-h+sAhif) 



/ dtAh{t)Cov^^{(T{t)J{a)) (3.6) 
Jo 
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and then 

Cov^^{a{t),f{a)) = ^Jf,n(.da)f,f,{da'){a{t)-a'{t)){f{a)-f{a')) 
= 2 / M(da)M(da')(/('^)-/('^')) 

J(T(t)=+,(T'{t) = - 

= = +) lih {<T{t) = -) [llh {fWit) = +)-llh {fWit) = -)] 

Conditionally on <T(t), the spin on [0, t] and are independent. So if we take the product of 
two couplings as in Lemma 13.41 on these two time intervals, we obtain a monotone coupling ^ 
of /X (.|<T(t) = +) and /i (.|cr(t) = — ) with the property that 

= +,cr" = -) = n ^^t>h^{o■ = e\a{t) = e)iX(t,p),h,x{o- = = e) 

e=± 

= IJL{a = +\(T{t) = +)ijL{a = -\cT{t) = -) . 

Furthermore, 

M(/k(t)=+)-M(/Kt) = -) = ^{f{a+)-f{a-)) 

^ M/(^+ = +,^~ = -)(/(+) -/(-)) 

so we conclude that 

Cov^J^(t),/(tT)) ^ (a = +)/.^(^ = -)(/(+) -/(-)) (3.7) 

which, in combination with (j3.6p . proves the statement (j3.3p when vr = 0. It is not difficult to 
derive an explicit constant as ip{a = +) = ip{a = —) = exp (— /3A) /2, the density of /i with 
respect to (p being bounded. Now we address the case of periodic imaginary time boundary 
conditions. We call f/ff = Hh {-{cr^O) = e,cr{f3) = e') the measure /x^, conditioned on <t(0) = £ 
and <t(/3) = e' . It is clear that 

/^r = (3-8) 

where 

and /^(o") = l{a-(o)=(T(^)=e}) for e = ±. Note that and — /~ are increasing, and their 
amplitude is ef^ (+) — ef^ (— ) = 1. When we take derivatives, it is a consequence of (|3.3|) for 
free imaginary time boundary condition that 



dPh+sAh 




ds 


s=0 



-^fJ-h+sAhif^ 



M(r) 



s=0 



> C 



+ -T-fJ-h+sAhi-f ) 

ds 

dtAh{t) 



s=0 



-m2 



(3.9) 



Now we take derivatives in (j3.8p : 
d 



per 



^g^^h+sAh 



.(/) 





dph+sAh 


- s=0 


ds 



(/) - I^H- if)) 



s=0 



+P-^f^htsAhif) + (1 - P) ^f'h+sAhif)- 



(3.10) 
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As the imaginary time boundary condition can be realized by adding to h an additional field 
eAks, where ks = l[o,5]u[/3-<5,/3]> ™ the limit A — > +00 and 5 — )■ (in that order), the last 
two terms are non- negative according to (j3.3p for vr = 0. So in order to conclude the proof of 
(|3.3p for periodic imaginary time boundary condition it is enough to provide a lower bound on 
f^h^(f)-f^r(f)- We know that 

f^h if) = lim , lim Mh+eAfc, (/) • 

Furthermore, according to (j3.3p for vr = 0, for any 5 > 0, ^4 G R i-> fih+Aks (/) is a increasing 
function and, on the interval [—1/(5,1/(5] its derivative is not smaller than c'(5 (/ (+) — /(—))• 
This proves that 

hm hm fih^Aksif) ^ c (/(+)-/(-)) (3.11) 

A— s>+oo A— >-+oo 

for any 6 > 0, where c' does not depend on 6, and consequently //^"^ (/) — fij^~ (/) satisfies the 
same lower bound. Putting p. lip with ()3.9p into (|3.10p we obtain (|3.3p for periodic imaginary 
time boundary condition. □ 

Proof of Corollary \3.3i Fix e <C 1 and assume first that /, g are increasing functions of a single 
spin, measurable w.r.t. J^r+e, for some T € [0,(3], where J't is the ci-algebra generated by 
{a{s)}s<t- We prove that 

Cov if g\TT) ^ -Ce^ll/ILIblloo (3-12) 

where 

Cov if, g\J^T) = ;u (/^I-Ft) - /i (/I J-t) /i (sI-^t) . 

Indeed, call the distribution /i conditioned to {cr (t)}t^T and to the event that cr has at most 
one flip in [T, T + e]. The distribution is completely described by the law of the time of the flip 
(+00 if no flip). But / and g are both monotone functions of this random time (both increasing 
or decreasing), so it follows from the FKG inequality for distributions on the real line (Lemma 
16.2 in the lectures notes by Peres [8]) that u (fg) ^ z/ (/) u (g). But the total variation distance 
between 1/ and /x (.|cr (t) , t ^ T) is less than Ce^, which proves (13.12p . 

Now we take two arbitrary increasing functions f,g of a single spin and choose T = (3 — e. 
When we apply the standard formula for conditional covariance together with (|3.12p we get that 

Cov(/,5) = ii{Cow{f,g\FT)) + Co^{lJi{f\TT),li{g\TT)) 
^ -Ce^ ll/lloo II5II00 + Gov (^ (/I Jt) , ^i (sI-^t)) 



where /i (/jj-y) and ;U {g\FT) are increasing functions with infinite norm less than that of 
and II^IIqo, respectively. By applying (j3.12p repeatedly with T = (3 — ke, k = 2, .../3/e we 
conclude that 

Cov(/,5) > -Cf3e\\f\\^\\g\\^ 
and the claim follows by letting e — )• 0. □ 



4. Glauber dynamics for the Quantum Ising model 



4.1. Definition of the generator and the semi-group. Here we define the Glauber dynam- 
ics for finite graphs and establish some preliminary properties. 
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The dynamics consists in resampling spins locally, according to the field generated by their 
neighbors. Given the graph G and the parameters /3 > 0, A ^ 0, /i : F x [0, /3] i-> R integrable, 
we call /Li = lJ^G;i3,h,\ the Gibbs measure on G with corresponding parameters and 

/ig = f,{.\ay = Py,yyeV\{x}), VpGS^,VxGF. 

Note that fj,x takes into account both the field h and the boundary condition on V. According 
to the DLR equation, is the measure obtained by taking a = ponV\ {x} and cTx according 
to the Gibbs measure on {x} with field + Ylyr^x Py '^^'^ interpret ^x as a kernel, since 
for each p € TX and x € is a probability measure (furthermore, p i-^ fix is continuous, 

cf. (j3.2p ). Because fi§ is a conditional expectation, it is a contraction in (^). 
Next we define the generator of the Glauber dynamics by 

C = ^ifix-I) 

where / is the identity operator. This is clearly a bounded operator on {fJ'G;fS,h,x) and the 
associated Markov semi-group is 

P, = e'^ 

= 5:e-'l"l<^x(')" E ......... (4.1) 

n>0 ■ VI 1/ xi,...,x„£V 

According to its definition, Pt is a bounded operator on (/x). Equation (14. ip shows as well 
that, for any t ^ 0, is a convex combination of the iterates of Hx, and is therefore a Markov 
kernel that contracts L^. For any p G S^, we will write Pf for the probability measure on 
which, on bounded functions /, acts as Pf (/) = (Ptf) (p). If p is distributed according to a 
probability measure on S, we will write lyPt for the measure J Pfdv (p). 

4.2. Mixing time and spectral gap. In this section, in analogy with the classical situation, 
we prove some basic results that are useful to control the relaxation of the dynamics to the 
equilibrium Gibbs measure. 

Proposition 4.1. Let some G = {V,E) finite, (3 ^ 0, h : V x [0, (3] ^ H integrable and A ^ 0. 
Define 

T^ix = mi{t;^0:yp,rieJ:^,\\Pf -P^W^^^e-^}. 
Then Tmix < oo and, for any t ^ 0, 

sup IIP," - ^ sup WPf - P[>\\^^ ^ e-L*/^-'''J. 

Proof. The first inequality is a consequence of being invariant by Pt (in other words, is a 
convex combination of the P/'). The second inequality is classical consequence of 

dit) = sup ||Pf-P,^||^^ 

being sub-multiplicative, cf. |5|. It remains to prove that Tmix < oo, or equivalently that 
d{t) < 1 for some t > 0. This follows from (j4.ip once we remark that any fix^ ■ ■ ■ fix„ with 
{xi, . . . ,x„} = V gives a probability at least c"'^' to the uniform plus state, uniformly in the 
starting state p. □ 
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Proposition 4.2. Let some G = (V, E) finite, /3 ^ 0, h : V x [0, /3] ^ H integrable and A ^ 0. 
Define 

ir, ■ f Gov (/,-£/) 

gap(£) = mf — — — — r- — 

/GL2(M):Var(/)>0 Var 

where Gov and Var re/er, respectively, to the covariance and the variance under fi. Then 

(1) There exists c < oo depending on f3,X and \\h\\^ such that 

gap(/:) ^ c^l^l. 

(2) For any f e L"^ (fi), 

Var(Pt/) ^ e-2*g^P(^)Var(/). (4.2) 

Proof. The proof of the second point is standard. For the first one we refer the reader to the 
proof of Theorem 6.4 in the Saint Flour course [6]. □ 

4.3. Monotonicity. Now we address the question of the monotonicity of the dynamics. An im- 
mediate consequence of ()4.ip together with the monotonicity of the single site measure (Theorem 
13. 2p is the following fact: 



Proposition 4.3. Take p,T] G such that p ^ r], and h ^ h. Denote by Pt the semi-group 
corresponding to field h. Then, 

Pf ^ P^. 

stoch . 

According to the convergence towards the equilibrium measure, it follows that 
Corollary 4.4. increases stochastically with the field h. 

Remark 4.5. The same argument as above could be used to establish the existence of a grand 
coupling, but for this we would need to know the existence of a grand coupling for the family of 
single spin measures, given an arbitrary family of external fields. 

4.4. Censoring. For any A CV, we let 

= ilJ-x - I) ■ 

Now we consider a function A : R"*" — )■ V (V) with finitely many discontinuities at to = < ti < 
... < tn- We define the censored dynamics according to A by the kernel 

where k is the largest integer in {0, . . . ,n} such that t ^ t^, and Ai = A {tf). Of course, when 
A{t) = V for any t ^ we get Pa^i = Pt, the uncensored dynamics. The theory of censoring due 
to Peres and Winkler [8] also applies here. Remarkably their result on total variation extends 
also to variance and entropy. 

Proposition 4.6. Consider A, B : V iV) as above. Assume that, for any t ^ 0, vl (t) C 

B{t). Assume that v is absolutely continuous with respect to ix with du/dfj, G L'^{dfi) and 
increasing. Then, for any t ^ 0, both vPa-i ctnd vPB;t o.fe absolutely continuous with respect 
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to /i, their Radon- Nikodym derivative is increasing and vPB;t ^ T^PA\t- Moreover the following 
inequalities hold: 



Var 



Ent 



d {uPB;t] 
d{iyPB;t: 



dfj, 

WPB;t - ^IItv 



Var 



Ent 



d {l^PA;t) 
d{l^PA;t) 



^ WuPA-t-fJ-l 



TV 



(4.4) 

(4.5) 
(4.6) 



Remark 4.7. A special case satisfying the assumptions of the proposition is when is concen 
trated on the identical equal to + configuration. In that case we will write vPA;t 

i^Pb± 



P^ . and 



Bit- 



Following [8j we begin the proof with two lemmas. 

Lemma 4.8. Consider v some measure on TX , and assume that v is absolutely continuous with 
respect to fi with ^ being an increasing function. Then ufix is absolutely continuous with respect 
to and its Radon-Nikodym derivative is increasing as well. 

Proof. Let cr ^ r and assume that dv/dfi is increasing. We denote by <t* the spin configuration 
on y \ {x} equal to <t on y \ {x}, and by <t|. the spin configuration on V equal to ^ at x and 
to cr on y \ {x}. Then 

d{u^^x),_, _ d.^^.^^|^^^^|^.,dz.,_,, 



d/j, 



dfi 



/x(d^|a^)^(r|) 



dfi 



dfi 



(r) 



where we used, in the first line, the fact that the density does not depend on cr^, since that spin 
is resampled, and at the second line the assumption that dv/dfi is increasing, and finally the 
fact that the single spin marginal increases with the external field (Theorem 13. 2p . □ 

Lemma 4.9. Consider v some measure on TX , and assume that v is absolutely continuous with 
respect to fi with ^ being an increasing function. Then vfi^ -< v. 

Proof. Contrary to [8j the set of single spins configurations is not totally ordered. For this 
reason we use an alternative proof based on the FKG inequality for single spin measures. Let / 
increasing. We start with 

y{f)-{^l^.)U) = ^(/-/^x(/)) 

and condition on a* , the spin configuration outside x. We have 



l^x 







1 



f 



du 
d/j. 



J^x [ d^J, 



l^x 



if) 



where in the second line we use the FKG inequality for a single spin (Corollary 13.3 



□ 
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Proof. (Proposition 14. 6|) . Lemmas 14.81 and 14.91 together with formula (j4.1|) imply that vPa-^i 
and vPs^t are absolutely continuous with respect to /i, that their Radon-Nikodym derivative is 
increasing and also that vPb^i ^ vPA;t- Now we prove the inequalities. For the variance, we 
remark that 




d{uPA;t) d{vPB 



< Var('li^V%ar^^(^^-^*^ 



d/i d^ 

1/2 



which proves (|4.4p . For the entropy we recall that 

Ent(/) = sup{^(/5):5with/x(e9) = l} 



therefore 



^.U'jp^) = (.p.,) flog ^C^-) 



dfi J \ d^ dn 

dil^PA;ty 



dfi J ' \ d/i 

< Ent , 

dfi 

Finally we recall for completeness the proof of (14. 6p : 



5. ISING MODEL ON REGULAR TREES 



□ 



5.1. Notation and main results. In the following we specialize to the case where the under- 
lying graph is T^*, the rooted regular tree with 6^2 children at each node except the leaves, 
and depth / ^ 0. We always denote by r the root of the tree. When z E T^* we denote by \z\ 
the depth of z, that is the graph distance to the root. We say that z is a leaf if \z\ = I. Given 
T € and ^ is a subset of we define //^ as the Gibbs measure on T,^ with boundary 

condition r acting as an additional local field at z G ^ given by Yly^Ay^z"^ iv)^ where y ^ z 
means that y,z are neighbors. When r is identically equal to plus/minus we will simply write 
/i^. When A = T'I, we will write /uj". 

Definition 5.1. We say that the parameters /3, A ^ 0, /i € ([0)/?]) are inside the uniqueness 
region if the boundary condition on the leaves of Tj* does not change the marginal distribution 
at the root in the limit / +oo. Because of stochastic domination, this is equivalent to 

lim WfJ-tio-r e ■) - n7 {ar e ■)\L^ = 0. (5.1) 
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Next, following we define the exponents 7 and k as follows: 
Definition 5.2. Given /3,A ^ 0, /i G ([0,/?]) we let 

7 = supmax ||//^ (<T^ G •) - /x^ (cTj, G •)||^^ (5.2) 

where the maximum is taken over all subsets A C T^^, all vertices y on the external boundary 
of A, all boundary configurations T],fi ^ T^'^'^ that differ only at y and all neighbors z ^ A oi y. 

Given a boundary condition r G S""""", we let k (t) be the infimum of k ^ such that there 
exists ^0 ^ such that, for any regular subtree T C with root x and uniform depth, and for 
any z £ T with \z — x\ ^ Iq, 

/i^^ • Ikx = +) - • l^i: = -) ^ /^'^'^ (5-3) 

From Proposition 13.11 we know already that 7 < 1. It follows from a recursive coupling 
argument that /t ^ 7. In complete analogy with the classical case, we prove the following 
results. 

Theorem 5.3. (Decay of correlations) Let (3,X^ 0, h £ ([0,/3]). 

(1) K{ + )^l/b. 

(2) If l3,X,h are in the uniqueness region, then for arbitrary t G k (t) ^ 1/6. 

Theorem 5.4. (Fast mixing) Let /3, A ^ 0, ft. G L"*^ ([0, /?]). Fix t G and assume that k (t) 
is such that k (t) 76 < 1. Then the following holds. 

(1) The spectral gap of the Glauber dynamics on T\ with boundary condition r is greater 
than a positive constant which does not depend on I 0. 

(2) The mixing time corresponding to the above dynamics is at most CI where C < 00 does 
not depend on I ^ 0. 

5.2. Conditional spin distributions. Working with trees leads to major simplifications in 
the structure of the Gibbs measure. In particular, given any 2: G T, the restriction of a Gibbs 
measure onto the subtrees of z, given cr^, is a product measure. As a consequence, n is fully 
characterized by 

(1) the marginal distribution fi [a^ G .) of the spin at the root r, 

(2) and the conditional marginal distributions /i {a^ G .\a^-), z G T\{r}, where z~ denotes 
the ancestor of z. 

Remark 5.5. The marginal distribution /i {a^ G .) can be viewed as a conditional spin distribution 
if we add a ghost ancestor to the root r with constant spin ctq = 0. 

For latter purposes we need that (conditional) measures form a vector space. So we introduce 
the set M of finite signed measures on (Ti,B(Ti)), where B(Ti) is the Borel cr-algebra associated 
to the Skorohod topology on S. Signed measures have a Hahn- Jordan decomposition into their 
positive and negative parts. This means that, for any G Al, there are two disjoint sets P and 
N (unique up to /x-negligible sets) such that ^ gives non-negative weight to every Borel subset of 
P, and non-positive weight to every Borel subset of A^. Let /i"*" = ^(. H P) and /i~ = — /i(. fl A^). 
Then fi = fi^ — fj.~ . We also consider \fi\ = ^'^ + ^~ , a positive measure, and recall that 

IImIItv = = \ sup ^^{f ) = I sup (/X (A) - /X (^^)) (5.4) 

defines a norm on Ai, that turns A4 into a Banach space as S endowed with the Skorohod 
topology is a Polish space. 
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For convenience, we call Mq the subset of A4 made of all finite signed measures p with 
= 0. Both M. and A^o are vector spaces. We also call M.+ the set of finite positive 
measures and the set of probability measures on S. 

Finally, we call X the set of functions from S to A^, Xq the set of functions from S to 
A4o, and similarly X+^i for the set of functions from S to In particular, any marginal 

(resp. conditional marginal) distribution belongs to A^+,i (resp. ^^+,1), and any difference of 
marginal (resp. conditional marginal) distributions belongs to A4q (resp. ^o)- For notation 
consistance, for any p & X and any cr G S, we denote by p"^ {■) the corresponding signed 
measure on S. We will consider later on the problem of defining a norm on X (see (j5.14p and 



5.3. The resampling operator and the cavity equation. Given r] G ([0, /3]), p ^ M. and 

pi, . . . , G we define Rpi,...,p^{p) G 7W by 



We give an interpretation to Rpi,...^p^{p) in an important special case. 

(1) Consider some point of the tree T with 6 + 1 neighbors. Conventionally we denote its 
spin by ctq. We call rj the spin of its ancestor and cri, . . . , cr^ the spins of its children. 

(2) Sample ctq from p, the conditional distribution of ctq given the neighbor spin rj. 

(3) Sample each cr j independently, according to pf , the conditional distribution of cr^ given 
the parent spin ctq. 

(4) Sample again ctq according to the single site distribution with field /i + 17 + <ti + • • • + <T{,. 

If p and the pi correspond to the conditional marginal distributions of the Gibbs measure with 
field h, then p is stable under the resampling operator according to the DLR equation. In other 
words, p = Rpi,...,p^{p) when p is the conditional distribution of ctq given that the spin above 
is rj, and the pf are the conditional distributions of cTj given ctq. More generally we have: 

Definition 5.6. Given G X, r] G L^([0, /3]), we say that p G Ai satisfies the cavity 

equation with parameters pi, . . . , pb and t] if 



5.3.1. The resampling operator is a contraction. From the definition (15. 5p it is clear that the 
applications /? G H' Rpj^,...,p^{p) G A4 and pi G X 1-^ Rpi,...,pi,{p) £ M are linear. Now we 
consider their operator norm. 

Proposition 5.7. Let 7 be the constant in Proposition I3.il corresponding to the value of M 
given by M = max (A, |j^||i + /3{b + 1)). Then, 

(1) For any pi, . . . ,pi, G ^1,+, any 77 G S, the application i?pi,...,pj is a contraction on M., 
and a ^-contraction (i.e. its operator norm is at most j) when restricted to A4o. 

(2) For any pi,. . . ,pi, G -^1,+, any r] G T,, the application I — Rp^^...^p^ is invertible on A4o 
and its inverse 



dsn])). 




(5.5) 



(5.6) 




(5.7) 



has operator norm not larger than (1 
(3) For any p e M, pi, . . . , pb e X , 



7) 



-1 




TV 




(5.8) 
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Proof. (1) We write p = p'^ — p the Hahn- Jordan decomposition of p into positive, mutually 
singular measures. We have 

||-^pi,...,Pb (p) llxv ^ ll^piviPb (p ) IItv ~'~ ll'^pi'-'-.Pb (p )IItv 
But (p^) is obviously a positive measure, and its mass is 

Rl_^^ (p+) (S) = j p+(d^o)pr'(dTi) • • • pV'{^<t,) 

therefore ||i?pi,...,pj, (/9)||rpy ^ IIpIItv • ^^"^ assume that p G A^o, that is, = 
The same calculation as above shows that Rpi,...,p^ (p) € A^o- Now, given 
(To,?7 S S we consider the probability measure 

= j pT{^^l)---p?{^<^b)t^P,h+r,+.,+-+.,^. 

It is immediate from (j3.ip that 
Consequently, 

Rl,...,p, (p) i^) - Rl,...,p, ip) (An 

_ f p+{dao)p-{dcT',) 

^ 27||p||tv 

since p'^{'^) = = ||p||tv- This shows that Rp^,...^p^ is a 7-contraction on ^Ao. 

(2) It is an immediate consequence of the first point that p G ^Ao 1— )• p — Rpi,...^p^ (p) € Mq is 
invertible, with the given inverse. The computation of the operator norm of the inverse 
is immediate. 

(3) This follows at once from the definition of total variation distance : 

|1;u||tv = sup p{f)/2. 

□ 

5.3.2. Monotonicity of the resampling operator. In this paragraph we examine the monotonicity 
properties of the resampling operator. We recall that, when /i,z^ G 7Wi^+, we say that fi is 
stochastically larger than when, for all / : S ^ R increasing, n{f) ^ t^if)- We generalize 
this notion by saying that p G A^o is stochastically positive when p{f) ^ 0, for all / increasing, 
which we write p y 0. Note that this has nothing to do with "p is a positive measure", which 
itself is the same as p G 

Remark 5.8. Let p, v two probability measures and assume that p is stochastically larger than v. 
Let p = p — v and call p^ the positive and negative parts of p. We claim that a coupling for the 
positive and negative parts of p yields easily a monotone coupling of p and i^, which is optimal in 
the sense that it realizes total variation distance. In other words, there always exists an optimal 
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and monotone coupling between two ordered probability measures. An explicit construction of 
this coupling is given below. By definition of the total variation distance, 

/9+(s) = = ||/i-i/||Tv- 

Assume that v and consider ip a monotone coupling of the stochastically ordered probability 
measures — i/||^y/5^. Note that — p^) = {y — p~) is a positive measure with mass 1 — — 
i/||tv- We call ip the law of {a, a) when a ^ {fi — p^)/(l — \\p — i^Htv)- Then, 

IIm - z^IItv93 + (1 - \\p - i^IItv) 

is a monotone and optimal coupling for (fJ.,!^). 

We begin by observing that Rpi^...^p^ {p) is increasing in all parameters p G ■M+,i, pi, ■ ■ ■ , Pb € 
<Y-|-i and T] such that for all i G {!,..., 6}, pf° increases stochastically with the boundary 

condition ctq € S. More precisely, -Rfi,...,pj (p) >- Rpu...,Pb {p) if P G ; pi, • • • , Pb, Pi, • • • , G 

Af+ 1 and r],i) G ([0,/3]) satisfy fj ^ r] and pf° >- pf°, i = 1 ... 6, for any (Tq G S. Also 
Rpi,...,Pb (p) >- Rpi,-,Pt ip) P^P^ -^1,+) Pi, • • • ,P6 e and T7 G ([0,/?]) are such that p >~ 
p and p""" increases stochastically with ctq G S, i = 1 . . . 6. Both statements follow immediately 
from equation ()5.5p and Theorem 13.21 

It appends that one can compare the stochastic increase in the parameter pi with the norm 
of the difference: 

Proposition 5.9. Let ij ^ Ti. Let c > and T be the constants in Proposition \3.1\ and Theo- 
rem \3.S\ corresponding to /? and M = max(|j/i||i + + 1),A). Let / : S — )■ R 6e increasing. 
Assume that p G A4+^i, pi G Xq with p^" >- for all ctq G S and let p2, ■ ■ ■ , Pb G ^+,1- T/ien 

(1) 

\\RI,-,pM\ty ^ r J p{dao)pricT*l). 

(2) 

Rl,...,pMif) ^ c{f{+)-f{-)) j p{da,)p1^{a.l). 

Proof. (1) The assumption that p^" is stochastically positive means exactly that its positive 
part is stochastically larger than its negative part. We consider therefore ^"'^ a monotone 
coupling of these two parts: 

||^pi,...,P6(p)|Itv 

= \ I ' '^r)Pr(d'T2) • • • pVidab) [Pn+r,^,^ if) - Ph+r,+.- if) 

^rj p{dao)d^'''{cT+,CT^)\\a+ -a^Wi 
= tJ p{dao)d^''''{a+,a^){at -a^)*l 

= T j p(dao)pr('^«l) 

(2) Let F{(Ti) = J p2°{dcT2) ■ ■ ■ p^°{dab)ph+r,+cTi+-+(Ttif), then ai H> F (ai) is non- 
decreasing and moreover, for every ^ af, 

F{at)-F{cT^) ^ c(/(+)-/(-))K-^r)-l 
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according to (|3.3p in Theorem 13. 2i When we introduce the same coupling of the positive 
and negative parts of p^'' as in the former paragraph, we obtain the sought lower bound. 

□ 

5.3.3. The cavity equation has a unique solution. 

Theorem 5.10. Let pi,...,pi, G Af+^i and r] £ L^([0, /?]). The cavity equation i5. 6\) has a 
unique solution in M^^i, that we call <l>pi,...,p5. Furthermore, this solution satisfies: 

'^l,P.,...,Pb - '^l,...,Pb = (^-^^i,P.,...,pJ''«^i-Pi.P.,...,P.(^^i,...,pJ^ Vp-iGAfi,+. (5.9) 

Moreover, the solution to the cavity equation $pi,...,pj increases stochastically with its parameters 
pi,...,Pb G -^+,1 increasing with the boundary condition, and r] € L^([0,/3]). More precisely, 
>- ^Pi,-,P6 -if Pi,- ■ ■ ,Pb,Pi,- ■ ■ ,Pb G X+^i and T],fi € ^^([O,^]) satisfy fj ^ r] and 
pf >- p^° , i = 1 . . .b, for any ctq € S, with both and being stochastically increasing in 
ctq G S. 

Proof. We consider first the case when pi £ are constant, i.e. pf does not depend on cr G E, 
for alH = 1 . . . h. Then, i?pi,...,p^ (p) depends uniquely on the mass of p G A4. Consequently the 
unique solution to the cavity equation in A^+,i, for this parameters, is Rp^,...,p^ (u) where is 
an arbitrary probability measure on S. Now we consider pi,. . . ,pj, G such that the cavity 
equation has a unique solution <I>pi,...,pj together with an arbitrary pi G and p G A^+,i. We 
use the shorter notations R = R%,...^p^^, R = R^i^,,„p^ and ^ = ^p^,...,p^. We have 

p = R{p) ^ p- <^ = R{p) - R{<i>) 

^ p-$ = ii(p-$) + i?($)-i?($) 
^ {I - R) {p-<^) = {R- R) ($) . 

But p — <I> G A^O) where / — ^ is invertible because of Proposition 15.71 Consequently there is 
a unique solution to the cavity equation with parameters pi, p2, ■ ■ ■ , Pb G ^^+,1 and it is given 
by formula (j5.9p . In the same way, when we change any other of the pi G <^+,i the solution to 
the cavity equation still exists and remains unique, so we can change all pi,. . . ,pi, to arbitrary 
elements of X+^i- 

The monotonicity of the solution to the cavity equation follows at once from the fact that, 
for any probability measure v G A^i,+, 

n....Pb = ir(^p"i....,pj'H (5-10) 

(according to the first point of Proposition ISTfl) together with the monotonicity of Rpj^,...,pi, (p) in 
all parameters p G pi, . . . , pfe G and r] such that for alH G {1, . . . , b}, pf° increases 

stochastically with the boundary condition ctq G S. □ 

5.3.4. Existence of a fixed point for the cavity equation. We consider 

= P^^^icTrG.) (5.11) 

for each n ^ and G S, where is the quantum Gibbs measure with parameters /3, h, X on 
the tree T„, with plus boundary condition on the leaves and field rj at the root. Observe that 
uq = 6^ and i/„,_|_i = ^u„,...,u„, the solution to the cavity equation with parameters Un, . . . , f„. 

Lemma 5.11. Let /3,A,M > 0. There is c> such that, for all h,r],f) with L^-norm at most 
M and all n G N*, 
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(1) ifv^v, 

{u^-i^:i){a.l) ^ c{v-rj).l. (5.12) 

(2) 

c ^ ^ ^ (5.13) 

dtp 

Proof. The first point is a consequence of Theorem 13.21 together with the DLR equation while 
the second one follows from the definition of fj,!^^ together with the DLR equation. □ 

Proposition 5.12. For every ij € ([0,/3]), the stochastically decreasing sequence of probability 
distributions Vn converge to a probability distribution vj^ on S. The conditional distribution 
Voo £ '^1,+ is a fixed point of the cavity equation, i.e. Voo = ^uoo,--;'^oc- 

Proof. The fact that the sequence is stochastically decreasing is an inductive consequence 
of the monotonicity of the solution to the cavity equation as z^i ^ i^o = <J+- Now consider the 
vr-system 

n = {{^(t) = +,Vt G /},/ C [0,/3] finite}. 

Obviously a (11) is the whole c-algebra corresponding to the Skorohod topology. Furthermore, 
any j4 G 11 is an increasing event and consequently Un (A) has a decreasing limit. On the 

1 77 

other hand, there is C < 00 such that the Radon-Nikodym derivative is uniformly bounded 
by C, for any n ^ 0, and consequently (fn) is tight. Indeed, if is the compact set of 
spin configurations with at most k flips, then Un (5^^) ^ C'v(^fc)- This proves existence and 
uniqueness of the limit v^d- The conclusion that z^oo is a fixed point of the cavity equation is an 
obvious consequence of the continuity of the solution to the cavity equation along its parameters, 
see (|5.9p and (|5.8p . Note that if we are not in the uniqueness regime, the cavity equation has 
another fixed point corresponding to minus boundary condition. □ 

5.3.5. Derivative of the solution of the cavity equation at the fixed point. Two natural norms on 
X are 



X = supIIp^'IItv (5.14) 
I d<^(r7)||p''||TV. (5.15) 



\\P\\i,x 

The first norm makes oi X a Banach space as M itself is a Banach space. As far as the 
second norm is concerned, we observe that when is absolutely continuous with respect to ^p, 

y|i,;, = i/dv.(r7)d^(a)|^|. 

Now we introduce the derivative of the solution of the cavity equation $pi,...,p5 along the first 
variable pi in the direction p Xq, at the fixed point pi = ■ ■ ■ = Pb = ^oo- The formula follows 
at once from (|5.9p : 

D^{p) = (/-iiL,...,..)"'^X,...,.oo(^^).Vpe^o. (5.16) 

Note that, for ah 77 G S and p ^ Xq, D'^ (p) € Mq. Consequently D {p) : t) ^ D'^ [p) G Xq 
and D can be seen as a linear operator on Xq. Our main theorem gives a bound on the norm 
of the operator norm of the fe-th iterate of D, denoted by , from the space {Xq, ||.||i,a') to 
("^Oi IMIoo.A'): 
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Theorem 5.13. There is C < oo that depends only on l3,h,X,b such that, for all p ^ Xq and 

for all k ^ 1, 

\\D'{p)\\oo,x ^ ^llplk^t. (5.17) 
In particular, the spectral radius of D is at most 1/6. 

Remark 5.14. The Krein-Rutman theorem [3] states that the spectral radius of a strictly positive 
and compact operator is an eigenvalue corresponding to a positive eigenvector. Although we do 
not use this theorem, it is remarkable that in the proof below the asymptotic direction of the 
stochastically positive conditional measure Vn — t^oo helps to control the norm of D. 

Proof. We first establish two preliminary results. Let p (z Xq and assume that p"'° is stochasti- 
cally positive, for every ctq € S. Then (i) in Proposition 15.91 together with (ii) in Proposition 
15.71 implv that, for any r] gT,, 



\D^{p)\\Ty ^ ^y"^2.(d'To)p'^«('^-l) 



^ V'(dcro)p"''(T.l). (5.18) 

where in the second inequality we have used (j5.13p . On the other hand, with the same assumption 
on p, if we take / : S ^ R increasing with /(+) = — /(— ) = 1, then 

D^{p){f) ^ cj ^{daQ)p^^^{a.l) (5.19) 

for some c > depending only on /3, h, A, b. Indeed, each term in the expansion of (/ — Ruoo,--;i^oo) 
leaves the set of stochastically positive p (z Xq invariant. So 

^ ^cj u^{dao)p''"{rT.l) 

where at the second line we use (ii) in Proposition 15.91 Inequality (I5.19P follows then from 
(f5?T3]l . Note that already inequalities (ISTTSj) and (jSTTg]) show that \\D{p)\\oo,x ^ C\\D{p)\\i^x 
when /o'^" is stochastically positive, for every ctq € S. 

The derivative D is closely related to the way converges to its limit z^oo- As a consequence 
of (]5.9p in Theorem 15. 10| 

n+l '^OO ^Vn,...,Vn l^oo ,■ ■ -,1^00 

^U„,...,Un ^!/„,...,!^„,!^oo ~ 

_ + 



-1 



1 ^00 VI ^00 ^00 VI ^00 



(I-R^ R'^ f^'' ) 



(5.20) 
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with IIt^IItv ^ C\\vn — VooW^ X f*^^ some finite C, since Rpi,....p,,{p) is multilinear and bounded, 
cf. (|5.8p . It follows by induction that, for every k ^ 1, 

Vn+k -Vco = b^D^{vn " i^oo) + rk,n as n ^ oo, (5-21) 

where ||rfc,„||oo,A' ^ C'fe||f„ - z^oo|lL,A'- 

Now we consider p ^ Xq such that /9'^° >- for every ctq G S. Without loss of generality we 
assume that ^t" ~ ^- F^^om (j5.20p we know that Vn — Voo is close to the image by bD"^ of 
y^-i ~ ^Sd as n ^ oo. According to ()5.18p and ()5.19p this implies the existence of c > such 
that, for all n large enough, for every / increasing with / (±1) = ±1, 

- v^){f) ^ csnY>\\v^-u'^\\n^y = c\\un-Voc\\oo,x- 



From ()5.18p and the assumption that \\p\\i = 1 it follows that 

CV 

sup||Z)^(p)||tv ^ C' = ^ 

According to the last two displays, for every n large enough, for all r/ € S, 

v v 



Call c' = c/ (2C"). As D preserves stochastic positivity. 



D 



k-2 



c'D{p) 



1 — l^oo - rk-2,n 

y-2 



c'D^-^ {p) 



is also stochastically positive {rk-2,n was defined at (I5.20|) ). It follows that, for all G S, 

1 ||j^?i+A:-2 — Z^oo|Ioo,A' + lkfe-2,n|loo ;f 



(d'-'-Y {p){a.l) ^ - 



Jjk-2 



for any n large enough. We can obviously find a subsequence of n along which ||f„_|_fc_2 — '^oolloo A" 
llt'n — i^oolloo A'' so taking liminf„ shows that 



D 



k-l 



According to (|5.18p we have 
D' (p) 

and therefore, for some different C < oo 



^ sup [ ^(dao) (d''-'Y\p) {a . 1) 



D'ip) 



oo.X 



(5.22) 



for all k ^ 1, all /5 E with = 1, such that p"''^ >- for every ctq G S. Finally we extend 

(I5.22P to any p ^ Xq with \\p\\i x ~ ^- Consider 



pV 
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Note that both p^" and p'^' are stochastically positive for any ctq G E. For instance, if / : E — )• R 
is increasing, then p^" (/) = (/ (+) — / (— )) ||/o'^"||xv + /''^" (/) ^ 0- According to the triangular 
inequality (note that ^ 2) and to (|5.22|) . 

and we are done. □ 



6. Proof of the main Theorems 



6.1. Proof of Theorem 15. 3L Without loss of generality, we only consider here the case x = 
in the definition of k. So we fix some I ^ and z G T^*. We denote by A; = \z\ the depth of z. 
Now we consider the subtree of T^* issued from r that contains z, and root it at z. This is a 
regular tree with 6 + 1 children at z and h children otherwise (except on the leaves) . Its depth is 
not uniform and ranges between min (/ — A;, A;) and / + k. The boundary condition is uniformly 
plus, except on the leaf r where it is + or — . 

Let (zo = r, zi, . . . , Zfc = -z) be a path from r to z. For any i G {1, . . . ,k — 1}, we call v- the 
conditional distribution fi (cr^. G 'If^i+i = ^, = i) and pf^ = P (f^yj G -Ic^.^^ = ^) where 
yj is the j-th children of Zj+i for j = 2, . . . , 6, if we consider Zi as its first children. We have. 



$ ± , i ^ k 

!^i_l.P(i-l,2) v,P(i-l,6) 



1 



and also, if we call f^'*^ the marginal distribution fi (cr^ G -IcTr = ±) ana pj = /x 
where yj is the j-th children of z for j = 2, . . . , 6 + 1, apart from zi, then 



±) and 







(note the 6+1 parameters in the cavity equation) . 

Now we claim that, given a neighborhood V C ^^1,+ of i^oo in the 
that min — i) is large enough depending on V, every P[ij) and lies in V. This is clear for 
P(ij) and as it is the marginal distribution at the root of a large subtree with identical plus 
boundary condition. As — ^ 7 ~ 

()5.2p . the same holds for u^^ . 

According to (]5.9p . for all e > 0, for all i with min(z,/ 
not on k, I, 



^ norm, for all i such 



^ where 7 < 1 as been defined in 
- i) large enough depending on e but 



D 



-)\\ 



oo,X 



loo, AT 



where D, the derivative of the solution of the cavity equation, at the fixed point along 
the first parameter, has been defined at ()5.16p . This clearly extends as follows. Fix j G N* 
and e > 0. Then, there is o = a(e,j) that does not depend on k^l such that, for all i with 
min (i, I — i) ^ a, then 



oo,X 



oo,X 



According to Theorem 15. 131 we conclude that, under the same conditions, 

C 



1+] 



i+i 



oo,X 



+ e 



loo, A" 
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where C < oo depends only on B,h,X. From lli^ii — u-,. 11 ^ ^ 7 ll'^,-^ ~ II v it follows 

1 1 t-f- 1 t T" -L 1 1 oo , 1 1 t t 1 1 OO , 

that 



-,o -,o 



^ ( 17+^^ 
TV V 



which proves the first statement of the Theorem. 

Now we consider /3, h, A in the uniqueness regime. Let r E and denote nj the Gibbs 

measure on T^* with r acting as a boundary condition on the leaves. Similarly to the proof of 
the first point of Proposition 15.91 we have, for any y, z € T^* with y the ancestor of z, and every 



|/i;+ (cT^ e -IcTj^ = r]) - fij (cr^ € -IcTj^ = r/)| 



TV 



r/)| 



TV ■ 



X child of z 



This proves that, as before, every P[ij) and lies in a given neighborhood F C Afi^+ of i^oo in 
the _:t> norm, for all i such that min (i, / — i) is large enough depending on V but not on r. 



The rest of the argument is identical. 



6.2. Proof of Theorem 15.41 Let r E be such that (r) < 1. We recall that, from 

Theorem 15.31 the plus boundary condition always satisfies this requirement, and that, if /3, h, A 
are in the uniqueness region, we could choose any r. 



6.2.1. Spectral gap. Following [7j (see Section 3 and the proof of Theorem 4.3), for the proof of 
the first part of Theorem 15.41 it is enough to verify (see Equation (14) in [7]) that for all e > 0, 
for all /c ^ 1 large enough depending on e, for all / ^ k, for any g ^ L? (S, dip) and all 77 € S, 

Var^^(w(5('T.)|{^.}|,|=fc)) ^ ((l + e)K76)'/'Var^,(5(^.)) (6.1) 

where yii is the Gibbs measure on T^* with boundary condition r on the leaves and extra field r] 
acting on the root. In what follows, all the bounds will be uniform in 77. 

Let us denote K{crr, .) the Radon-Nikodym derivative of the measure fn (^fn (^cr'^. € ■{ {cTz}^^^^/.^ 

with respect to the marginal distribution of the root = fii{crr G •). The kernel K is constructed 
as follows: given cTr sample the spins at distance k from the root. Then take these spins and 77 
acting on top of the root as a boundary condition on and sample again cr^. It is useful to 
remark that K is uniformly bounded. The Cauchy-Schwarz inequality implies that 



s: Var^,, (m (g (cr^) | {cr^}, ^ Var^,, {g (cr^))^/^ . 
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Then 



1 



so (|6.ip would follow from 

\K{cTlcTl)-K{alal)f] ^ 2{{l + e)K^h)' 

From the boundedness of K we conclude that, for some constant C, 

.1 _3n T^r2 _3^^2" 



^ 4C sup \\K{al,.)i^ - K{cr'^,.)i^\\ 

^ 8Csup \\K{+, .)v - K{p, O^^IItv 
pes 



TV 



The DLR property, together with the definition of T at (|3.2p . imply that for ordered marginals 
the TV distance is comparable with difference of expectation of o" • 1. Consequently, 



^ C 



[\\Ki+,.)l.-Kip,.M^y 

-. Now, we argue that the former difference is bounded 



which is clearly maximum if we take p 
by 



z:\z\=k 



cr -, • 1 (Tr 



+ ) - (CT-^ • \ \cTr 



X 2r7' 



(6.2) 



For proving this we need only a slight adaptation of the proof of (ii) in Claim 4.4 in [7]. Consider 
two spin configurations that differ at a single position z with j^l = k. We can easily 

construct a coupling of /i/ [cTt G 'I {c^ = ^2}|;;|=fc) ^'^d /i/ ^cr^ € - j {<t^ = ^^jj^i^^^ for which 

the two variables differ with probability at most '^^~^V ||^' — ^||]^ according to the definition of P 
in Proposition 13. II and to that of 7 at (|5.2p . Now we consider a monotone coupling of {cr^}|^|^^ 
corresponding to the conditions cr,. = it. By applying the former coupling to an interpolating 
sequence between these spin configurations as in proof of (ii) in Claim 4.4 in [7], we conclude 
the proof of ([Ol . 

Finally, /i/ (a^ • Ijcr^ = +) — ^/ (cr^ • Ijcr^ = — ) ^ (\ ^ ef^ for all A; ^ large enough 
according to the definition of k. The proof of the first part of Theorem 15.41 is complete. 
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6.2.2. Mixing time. We first establish an intermediate step. We fix I and consider the Gibbs 
measure on with boundary condition r on the leaves of T^" and additional field 77 € S at 
the root. Now we turn to the dynamics. Consider as a starting configuration the identically 
plus configuration. Similarly we could consider the identically minus configuration. Let Pj^ 
be the distribution of the spin configuration at time t for the dynamics. We denote by hr {t) 
the Radon-Nikodym derivative of with respect to fi. Note that hr (0) = !{+}///(+) and 
hr (t) = Pthr (0) because Pt is self-adjoint in (^u). 

We also define Tr = sup^minjt : Var(/i7 (t)) ^ 1}. Similarly, given x G T^" \ {r}, we consider 
the dynamics censored everywhere except on the subtree rooted at x and call h^ (t) the resulting 
Radon-Nikodym derivative with respect to the Gibbs measure conditioned on being plus outside 
the subtree of x. We then define = min {t : Var (t)) ^ 1}. 

Now we prove that there exists a constant to independent of / such that 

Tx ^ max Ty + to. (6.3) 

y child of x 

For this purpose we use censoring (Proposition 14. 6p together with the first point of Theorem 
15.41 For simplicity we only consider the case x = r. We censor for time t = maxy child of r Ty the 
root and then run the uncensored dynamics for an extra time to to be determined later. Let 
hr (t) be the Radon-Nikodym derivative of the corresponding distribution at time t. Then, 

Var (/i;? (t + to)) ^ Var (/i^' (t + to)) 
= Var (Pto/i;? (t)) 
^ e-2s^P>^*0Var (/ij'(t)) . 

Now we prove that Var (hr (t)) is bounded uniformly in /. By construction we have 

^ ' y child of r 

where a^j^^ is the restriction of a to the subtree rooted at y. Finally, we remark that 

v2^ 



\sx{h'f{t)) ^ I^[{h7{t)) 



fi{ar = +)fi{^{h';! {t) f\ar = + 



^('^'^ = child of r 

^ n [VarM.k.=+)(^+WK,))+l] 



uiar = +) 

' y child of r 

which is smaller than 2^/// (cr^ = +) according to the definition of t. If we take to such that 

2'' exp(-2gap xto) ^ ^ 

^{(Tr = +) ^ 

then to is bounded uniformly in / and we are done. In conclusion, we have shown that Tr ^ /to- 

Remark 6.1. We observe that the recursive inequality (16. 3p is analogous to the one obtained in 
[7] (see Lemma 5.8 there) for the logarithmic Sobolev constant. In our context, this constant is 
easily seen to be infinite because of configurations with an arbitrary large number of fiips. 
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We know have to consider the dynamics starting from an arbitrary spin configuration ^. 
Choose t = Tr + cl for some c > to be chosen later on. We have of course 

||P,^-;u||tV ^ ||P/-P+||TV + ||^'t+-/i||TV. 

On one hand, 

IIP+-/.IITV = 11/^7 (t)-l|lLi(^)^ll/i?W-l|lL2(^)^e-s-P(*-^^)^e-^^P-' (6.4) 

according to the definition of T^. 

For the remaining part \\P^ — Pj^|[tv we use a coupling argument. Call s = + cl/2. We 

consider ^'s a monotone coupling of Ps -< Pft and denote ((t,<t+) its variables. According to 
Markov's inequality, we have 

^'^(3x,cr^«l ^ cr+.l-e) ^ 1'^(cr^ • 1 ^ cr+ • 1 - e) 

X 

^ e-i^P+(^,.l)-P|(^,.l) 

X 

^ e-i^P+(^,.l)-P-(a,.l) 

X 

^ ^(||P+-^||tv + ||P7-/^IItv). 

Now we take e = h~'^K According to the last display and to (|6.4p the probability cr^; • 1 ^ 

cr+ • 1 — e) can be made arbitrary small by taking c large. So with high probability under 
^s, at not position the spin configurations differ by more than e in ([0,/3]) distance at time 
s = Tr + cl/2 = t — cl/2. In the remaining time cl/2, we update the spin configurations at the 
same positions according to an optimal coupling of the marginals. Using Proposition 13. 1 1 we see 
that, as long as the spin configurations differ at every x by at most e in ([0, /3]) distance, each 
update put locally the same spin with probability at least 1 — Te. Therefore we conclude that 
\\P^ — P^\\i[y can be made arbitrarily small if c is large enough. This concludes the proof that 
the mixing time is bounded by Cl. 
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